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The loss rate of linear momentum from a binary system composed of compact objects (radially 
falling towards each other under mutual gravitational influence) has been investigated using the 
multipolar post-Minkowskian approach. The 2.5PN accurate analytical formula for the linear mo- 
mentum flux is provided, in terms of the separation of the two objects, in harmonic coordinates, 
both for a finite and infinite initial separation. The 2.5PN formulas for the linear momentum flux 
- - - are finally used to estimate the recoil velocity accumulated during a premerger phase of the binary 

evolution. 
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I. INTRODUCTION 



Gravitational waves from coalescing binary systems carry away energy and angular momentum of the source. For 
asymmetric binaries (composed of objects of unequal masses and/or with nonzero spins), there will also be a net loss 
of the linear momentum from the source. As a consequence, the center-of-mass of the source will receive a recoil in 
the opposite direction. This recoil accumulates until the two objects of the binary merge to form a single object and 
the source stops losing linear momentum. At this juncture, the remnant of the coalesced binary moves with a non 
zero kick speed along a straight line path in space. For a more detailed discussion on the phenomenon of gravitational 
wave recoil, see Ref. The phenomenon of gravitational wave recoil is extremely important in various astrophysical 
contexts such as the formation and growth of super massive massive black holes at the centers of galaxies. If the 
^ . recoil velocity of the remnant of the coalesced binary is more than its escape velocity from the host, then the host 
^ will not be able to retain the remnant and models that grow the super massive black holes via successive mergers 
. from other black holes will not be favored 0- An accurate estimate for the recoil velocities associated with compact 
\^ ■ binary mergers can be used to address issues like observations of super massive black holes at the centers of most of 
^ ] the galaxies in the local universe Q or their apparent absence in globular clusters and dwarf galaxies or to predict 
the population of compact binary systems in globular clusters. 
I The importance of this phenomenon has been realized widely in as trop hysics community and there have been 
numerous analytical or semi-analytical j3-[l6| and numerical studies (T^-HJ to compute this effect. All these studies 
compute the recoil effects due to the loss of linear momentum from compact binary systems (which either have mass- 
asymmetry and/or have non zero spin) moving in quasi-Keplerian or in quasi-circular orbits. Numerical simulations 
for nonspinning black hole binaries moving in quasi-circular orbit [l"7l-l20j have shown that the recoil velocity can be 
5h ] of the order of few hundred kms""'^ while for spinning case [2T|-|2^ the recoil velocity estimates can reach up to few 
9^.- thousand km s~^. 

Although, head-on infall and the subsequent merger of two compact objects due to gravitational wave radiation 
reaction effects would be an insignificant astrophysical possibility, nevertheless it has been studied extensively using 
various analytical/numerical approaches. The motivation behind such a study is many-fold. To start with, due to 
the axial symmetry of the system, the two-dimensional problem of compact binary motion becomes one-dimensional 
and hence the treatment becomes simple. This also can act as a toy problem for comparing various analytical and 
numerical approaches in their most simplified versions. In addition to this, head-on collision can be considered as 
an approximation to the merger phase of the inspiralling compact binary evolution. Finally, as pointed out in [2^ , 
head-on collision studies can be used to remove the uncertainties in the direction of the recoil of the remnant. 

One of the earliest attempts to compute recoil effects due to the radial plunge of a test particle into a Schwarzschild 
black hole is due to Nakamura and Haugan [26| using the black hole perturbation theory. Using a close limit approx- 
imation method Andrade and Price [23 first computed the recoil effects due to head-on collision of two black holes. 
On the numerical relativity front, Anninos and Brandt (28| computed the recoil velocity due to head-on collision of 
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two unequal mass black holes. Some other (relatively recent) analytical/numerical works [lH, [2^, H^l compute the 
recoil effects taking in to account the asymmetry in mass and/or in the spin. As far as PN calculations are concerned, 
although, the recoil effects in a head-on collision case have not been investigated explicitly, one can use expressions for 
the linear momentum flux from nonspinning inspiralling compact binary systems moving in general orbits \4, la , jSj to 
write equivalent expressions for the head-on case by using the following transformations ^ (as suggested in [3ll |32|): 

X = zn, -V ~ in, r = z, v = r ~ z. (1-1) 

Here, z is the separation between the two objects (under radial infall) at a given instant and z is the first time derivative 
of z, giving the relative speed of objects at that instant. The most recent related PN work Q gives 2PN accurate 
expressions for the instantaneous part of the linear momentum and hence one can use the above transformations to 
write the 2PN expression for the instantaneous part of the linear momentum flux in terms of z and z. In the present 
work, we not only calculate the instantaneous part of the flux explicitly for the head on case to a higher order (2.5PN 
as compared to previous 2PN calculations) but also compute additional terms contributing at the 1.5PN order and 
2.5PN order (tail contribution) whose nature has been discussed in more detail in the next section. 

In the present work, we compute the 2.5PN accurate analytical expressions for the linear momentum flux, in 
harmonic coordinates, emitted during the radial infall of two nonspinning compact objects under mutual gravitational 
influence. We study the problem for two different situations based on the initial separation between the two objects. 
In the first case we assume that initially the objects are separated by some finite distance ( we call it case (a)) and in 
the other case we assume that the initial separation between them is infinite (we call it case (b)). Linear momentum 
flux as a function of the separation between the two objects at any instant of time for the two situations, case (a) 
and case (b), are given by Eq. (|4.1ip and (|4.12[) . respectively. We use these results to estimate the associated recoil 
velocity for the two situations. Since linear momentum flux expression (Eq. (|4.1ip ) involves some integrals (Eq. (|4.6|) ) 
which can only be evaluated numerically, it is not possible to give analytical PN expressions for the accumulated recoil 
velocity for case (a) and thus has been computed numerically. However, for case (b), a 2.5PN accurate expression for 
the recoil velocity is given by Eq. (|5.6p . A graphical representation of our results have been given in Figs. [l][2] We 
find that the recoil velocity is maximum for a binary with ly ~ 0.19 and is of the order of ~1.6kms^^ if we terminate 
our calculations when the two objects are 5 Gm/c^ apart. 

This paper is organized in the following manner. In Sec. [Ill we first write the general formula for the linear momen- 
tum flux in terms of the radiative multipole moments of an isolated post-Newtonian source. Next, we use relations 
connecting the radiative multipole moments to the source multipole moments, to express the linear momentum fiux in 
terms of the source multipole moments. Section Hill lists all the inputs that will be required for computing the 2.5PN 
accurate analytical expression for the linear momentum fiux. In Sec. IIVI we present the 2.5PN accurate analytical 
results for the linear momentum flux, in harmonic coordinates, for two situations (case (a) and case (b)). In Sec. IVl 
we show how the expressions for the linear momentum flux can be used to compute the associated recoil velocity 
accumulated till any epoch of the binary's evolution (within the validity of PN approximations). Finally, in Sec. I VII 
we summarize our findings and discuss the numerical estimates for the recoil velocity in the head-on case. 



II. THE POST-NEWTONIAN STRUCTURE FOR THE FLUX OF LINEAR MOMENTUM: 

CASE 



HEAD-ON 



The general formula for linear momentum fiux, in the far-zone of an isolated source, in terms of two sets of symmetric 
trace- free radiative multipole moments {Ul, Vl), is given in [ssj (see Eq. (4.20 ' ) there). The radiative moments, Ul{U) 
and Vl{U), are referred as mass-type and current-type radiative multipole moments, respectively, and arc functions 
of the retarded time U in radiative coordinates. Here, L = 11^2 ■ ■ ■ ii represents a multi-index comprised of / spatial 
indicies and U is given by U = T — R/c, where T and R denote time of observation and the distance to the source 
in radiative coordinates, respectively. At 2.5PN order, the expression for linear momentum fiux, in terms of radiative 
multipole moments {Ul, Vl), reads 
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In the above, denote the I''* time derivative of {Ul,Vl}, tijk denotes the Levi-Civita tensor with 

£123 = +1 and O (l/c^) indicates that corrections of the order 3PN and above have been neglected in the present 
analysis. The expression for linear momentum flux, in terms of radiative multipole moments {UL^VL), is not very 
useful unless we show how these moments are connected to the actual parameters of the source. Fortunately, the 
formalism for connecting radiative multipole moments to the source-rooted moments, with the PN accuracy desired in 
this work, has already been developed [SJ] using the multipolar post-Minkowskian approach [35l - [40| . In the multipolar 
post-Minkowskian formalism, Ul and Vl arc first written in terms of two sets of multipole moments. Ml and 5l, 
referred as mass-type and current-type canonical multipole moments, respectively. Next, these canonical multipole 
moments, Ml and Sl, are written in terms of six sets of multipole moments, II, JliWl, Xl^Yl, Z^^ referred as 
source multipole moments. The multipole moments, II and Jl, thoroughly describe the source and are referred as 
mass- type and current type source multipole moments. The other four, ^l, Yl and are referred as gauge 
moments as they do not play any role in a linearized theory and only become important at nonlinear level. Reference 
[11] exphcitly lists all the relations connecting {Ul, Vl) to {Ml, Sl) (see Eqs. (5.4)-(5.8) there) and those connecting 
{Ml, Sl) to {II, • • • , Zl) (see Eqs. (5.9)-(5.11) there). Using these relations one can explicitly write expressions for 
radiative multipole moments {Ul, Vl) (and hence the linear momentum flux at 2.5PN order given by Eq. (|2.ip ) in 
terms of source multipole moments {II ■ ■ ■ Zl)- Before we express radiative multipole moments in terms of source 
multipole moments, we would like to bring in to the notice the fact that, for head-on case current-type moments {Vl 
or Sl or Jl) would not contribute as they are proportional to the angular momentum, JT", which vanishes for the 
head-on case. This allows us to re- write Eq. (|2.ip . in a form specific to a head-on case, and it reads 
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It is evident from the above, that moments appearing at the lowest order in the PN series need to be known with the 
highest PN accuracy whereas those appearing at a higher PN order need to be known with smaller PN accuracy, e.g. 
in the present case we need Uij and Uijk to 2.5PN accuracy whereas Uijki and Uijkim need to be known with 1.5PN 
and Newtonian accuracy, respectively. Now, making use of Eqs. (5.4)-(5.7) and Eqs. (5.9)-(5.11) of [s^ and keeping 
in mind that current type moments vanish for the head-on case, we write C/l in terms of source multipole moments 
in a form specific to the head-on case, which read 
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In the above, angular brackets (()) surrounding indices denote symmetric trace- free projections. Here, M, is the 
total ADM mass of the source and rg is an arbitrary length scale and provides a scale for the logarithms in tail 
integrals. This length scale was first introduced in the multipolar post-Minkowskian formalism and enters the relation 
connecting the retarded time, U in radiative coordinate to the retarded time, u=t-r/c in harmonic coordinates, which 
reads 
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In addition, note the presence of two types of terms in above expressions: the first kind involves multipole moments 
at any given retarded time U and are referred as instantaneous terms and the other kind involves integrals over 
time, referred as hereditary terms that require the knowledge of multipole moments at any time U' = U — t before 
U . Further, the hereditary terms can be split into two parts: terms with and without logarithmic factors inside the 
integrals. Integrals with logarithmic factor are called tail integrals and those without logarithmic factor arc referred 
to as memory integrals. 

Since the linear momentum flux involves 1*^' time derivative of mass- type radiative multipole moments (Eq. (|2.2[) ). 



first we need to write U^"^ in terms of source multipole moments.^ In terms of source multipole moments, U^^ take 
the following form 
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It was argued and then shown in [33 [ (see Sec. II there for a detailed discussion) that the presence of in the tail 
integrals at 1.5PN order is due to our use of the radiative coordinates and will disappear if we insert U (given by 

Ea. (|2.4p ) back in expressions for Ul (same would be true for u'^j^^). Upon doing so we can write expressions for u'^^ 
in harmonic coordinates which now will be free from the arbitrary length scale, rp, and read 
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^ The memory integral is a time anti-derivative and thus becomes instantaneous when we take the time derivative of C/^. 
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Equation (|2.6p along with Eq. (|2.2p gives 2.5PN accurate expression for the linear momentum flux in terms of the 
source multipole moments in harmonic coordinates, in a form specific to the head-on case. Next, the resulting 
expression can be decomposed into two distinct pieces namely: the instantaneous contribution and the hereditary 
contribution whose nature has already been discussed above. The total linear momentum flux reads 
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where the instantaneous part is given by 
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and the hereditary contribution reads 
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Now, if we know how the source multipole moments are related to the actual source parameters, with PN accuracy 
desired in the present work, and we have a suitable machinery to compute the time derivatives of the source multipole 
moments, we can express the linear momentum flux in terms of actual source parameters. With this motivation we 
move to our next section where we shall provide all necessary inputs that will be needed for computing the 2.5PN 
linear momentum flux in terms of the source parameters. 
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III. INPUTS FOR COMPUTING THE LINEAR MOMENTUM FLUX: RADIAL INFALL OF TWO 

COMPACT OBJECTS 

As discussed in SecU in this paper we aim to study the loss rate of hnear momentum (through outgoing gravitational 
waves) during the radial infall of two compact objects under mutual gravitational influence. Unlike the case of 
inspiralling compact binaries in eccentric or circular orbits (where the motion takes place in a plane) , for the head-on 
case, the problem becomes one dimensional and thus the treatment becomes relatively simpler. For such sources, 
expressions connecting source multipole moments to the source parameters, with the PN accuracy desired in the 
present work, have been given in Ref. [32j.'^ Below we list all source multipole moments (in harmonic coordinates) 
needed for computing 2.5PN linear momentum flux in terms of the separation between the two objects at a given 
instant (z) and the first time derivative of z(z), giving the relative speed of objects at that instant (assuming the 
motion takes place along the z-axis).* The mass- type source multipole moments read 
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Here, rii is the component of the unit vector, h, along the direction of motion and 7 is our PN parameter and is related 
to the separation (z), between the objects at any instant of time, by 7 = (Gm/c^z). In addition to this, one would 
also need IPN accurate expression for mass monopole (while computing hereditary terms), which can be identified 
with the ADM mass (M) of the system and Newtonian order expressions for gauge moments such as the one related 
to monopolar moment W and dipolar moment Yi and are given as 
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Having expressed, the source multipole moments in terms of the parameters of the source, now we need to compute 
relevant time derivatives of the source multipole moments. With mass-type source multipole moments and other 
required moments given in terms of z and i, whenever a time-derivative is taken, terms involving z appear and thus 
one would need an expression for z in terms of z and z in order to write the linear momentum flux in terms of just z 
and i. Reference [13 lists somewhat general 3PN expression for z (in terms of z and z) which can be used to write 
related expressions in SH, MH and ADM coordinates by choosing appropriate values for the parameters, a and /3 (see 
Sec. IHA of [13 for details). However, for our present purpose we just need 2.5PN accurate expressions for the z in 



^ Reference [s^l provides a 2PN expression for the mass octupole moment (lijk) however for the present purpose we need it with 2.5PN 
accuracy and this additional 2.5PN correction is new to this paper (see Eq. Il3.1bl l). In addition, the moment, Yi, was not needed for 
the energy flux calculations at 3PN order but is needed here with Newtonian accuracy and is also new to this work (see Eq. Il3.2cl l). 

* Unlike Ref. |32|| . where expressions for energy flux are given in standard harmonic (SH), modified harmonic (MH) and Arnowitt, Daser, 
and Misner (ADM) coordinates, here we only make use of harmonic coordinates for all relevant formulas. However, in the appendix we 
show how one can obtain equivalent analytical expressions for the linear momentum flux and recoil velocity in ADM coordinates. 
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With, source multipole moments and z expressed in terms of z and z, we can compute all relevant time-derivatives 
of source multipole moments appearing in flux formula (Eq. (|2.8p - (|2.10p ) and then can use them to write the linear 
momentum flux (at least instantaneous part of flux since hereditary contribution shall involve computing the integrals) 
in terms of z and z. However, following (sil. [32j. we would like to write the expression for the linear momentum flux 
as a function of the separation of the two objects, alone. Also, we would like to compute the flux of linear momentum 
for two different situations: case (a) the two objects in the problem, initially separated by some finite distance, start 
falling radially from the rest, under mutual gravitational attraction, and case (b) a similar situation of radial infall 
but assumes infall from infinity. In order to write the linear momentum flux as a function of the separation of the two 
objects, we need an expression for i in terms of z, with a certain PN accuracy (here it should be 2.5PN accurate). In 
addition to this, z{z) is also sensitive to the initial conditions (case (a) and case (b)). At 3PN order, z{z) has been 
computed in [32| for the two different situations we want to explore in the present work and will not be reproduced 
here. We directly quote the result. In harmonic coordinates, 2.5PN expression for i, in case of infall from a finite 
initial separation (zi) is given as 



-V2cVl - s/7 



1 



jy 47 . 

+s" I - H \ V' 

' 8 2 32 



h - f 

2 4 



— V2VT 
15 



1 



9 



s 

2 4 



"^7^/2^. - 



27 



55 



7iy-\ 

8 32 



123 



59 . 

-ly ly' 

16 



(3.4) 



where, s = z/zi < 1.^ Related expression for the case of infall from infinity can be obtained by setting s = z/zi in 
the above and then taking the limit as — > c», and it reads 



1+7 



-2 + r)+r 



27 „ 



32 



15 



(3.5) 



With these inputs we now are in a position to write the instantaneous part of the linear momentum flux in terms of 
the separation between the two objects under radial infall. However, the computation of hereditary contribution shall 
require IPN expression for the trajectory of the problem.^ The IPN trajectory for the two situations (case (a) and 
case (b)) have been given in [11] (see Eq. (3.23)-(3.24) and Eq.(3.26) there) and we simply recall it here (with slight 
change in presentation). For case (a). 



3/2 



gis) 



1 Gm 



(3.6) 



where g{s) = /i(s) - /2(s), ho{s) = /i(s) + 9/2(5) and hi{s) = 9/i(s) + /2(s) with /i(s) = - s and /2(s) 

arcsin ^/s. For case (b) , the above expression reduces to 



V2z^/^ 



3 VGV?^ 



15 Gm 

T"?7 



1-^ 

2 



(3.7) 



We now have all the inputs to compute both the instantaneous and the hereditary contributions to the linear momen- 
tum flux, given by Eq. (|2.8p - (|2.10p . and have been computed in the following section. 



^ Note that at 2.5 PN order SH coordinates and MH coordinates arc equivalent. 

^ Note that, the 2.5PN expression for i has been obtained by adding Eq. (3.8) and Eq. (5.3) of [32l | (as was suggested there) and then 
truncating resulting expression at the 2.5PN order. 

Note that the leading order hereditary contribution occurs at 1.5PN order and thus computation of hereditary contribution at 2.5PN 
order shall only require IPN inputs. 
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IV. THE 2.5PN LINEAR MOMENTUM FLUX 



A. The Instantaneous Contribution 



Instantaneous part of the linear momentum flux, in terms of the source multipole moments and their time derivatives, 
is given by Eq. (|2.8p - (|2.9p . Expressions for the source multipole moments (Eq. p.ip - (|3.2p ') and the one for z (Eq. p.3p ). 
in terms of z and i, can be used to compute the relevant time-derivatives of source multipole moments algebraically 
as functions of z and z. Next, in order to express the source multipole moments and their relevant time-derivatives, 
solely as functions of z, we need to make use of expression for i given in Eq. p.4p - (|3.5p . depending upon the case 
we want to explore (case (a) or case (b)). Using, source multipole moments and their relevant time derivatives, 
solely expressed as functions of z, in Eq. (j2.8p - (|2.9[) . performing contraction of indices and truncating the resulting 
expression at 2.5PN order, we can write 2.5PN accurate expression for the linear momentum flux as a function of 
separation of the two objects {z). 



I. Case (a): Inf all from a finite distance 



The 2.5PN accurate expression for the linear momentum flux, for the situation which assumes the radial infall of 
two compact objects (initially separated by some finite distance z^), in terms of our post-Newtonian parameter 7, 
reads 
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(4.1) 



where 7 ~ (Gm/c^z) and s ~ z/ Zi < 1. In the above, note that the leading order contribution to the linear momentum 
flux is proportional to the parameter s and hence will vanishes for the case where initial separation is assumed to be 
infinite (zi — > 00 i.e. s 0). This is expected since the Newtonian order linear momentum flux is proportional to 
the 4*'' time-derivative of the octupole moment (lijk), which vanishes for the case of infall from infinity.^ However, 

for the case of infall from some finite separation the I^jl survives [Hj, and hence we see a finite Newtonian order 
contribution to the linear momentum flux. 



2. Case (h): Infall from infinity 



For the case of infall from infinity the related expression can be obtained by setting s = z/ Zi and then taking the 
limit as Zj — >■ 00 we obtain 
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(4.2) 



This was first noted and discussed in [26l | and can be verified easily. 
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B. The Hereditary Contribution 



The hereditary contribution to the hnear momentum flux, in terms of time-derivatives of the source multipole 
moments, is given by Eq. (|2.10p . Computing hereditary terms is relatively less easy as compared to computing 
instantaneous terms since it requires one to compute integrals over retarded time spanning over the entire dynamical 
history of the source. Now, since the leading order contribution to the linear momentum flux occurs at relative 1.5PN 
order we need to compute the hereditary effects only with relative IPN accuracy in order to achieve relative 2.5PN 
accuracy for the present purpose. Moreover, only first two terms of Eq. (|2.10|) need to be IPN accurate as the last two 
already contribute at 2.5PN order. In addition to this, in order to compute hereditary terms with accuracy desired in 
the present work, essentially we need to evaluate only three integrals, since integrals appearing in 2"'' and 3'"'* term 
of Eq. ()2.10|) are essentially the same. Below, we list the three integrals we need to evaluate (note r — >■ z) 

/(f(u-T), (4.3a) 

I^liu-rl (4.3b) 

I^%iu-r). (4.3c) 

As discussed above, Ii and I2 need to be IPN accurate whereas we need /a to be only Newtonian accurate. 



h = 



dT 



' 2z) 12 



h = 



dT 



, / CT\ 97 




1. Case (a): Inf all from a finite distance 
In this case, the integrals listed above can take the following form [s^ 



h = 


I dT 
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Inf — (u 
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h = 


I dT 

Ju(zi) 
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60 
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30 



I^'^'Ut) 



(4.4a) 
(4.4b) 
(4.4c) 



Note that for the infall from infinity case, when Zi — ^ 00, u{zi) = u{oo) = —00. With, required derivatives of the 
source multipole moments, expressed in terms of z, and the IPN trajectory (given by Eq. (|3.6p ) we can evaluate these 
integrals and they read 
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(4.5c) 



where, Intl(s), Int20(s), Int21(s),Int30(s), Iiit31(s),Iiit4(s), Int5(s),Int6(s), Int70(s), Iiit71(s), Int80(s), Iiit81(s),Int9(s), 
IntlO(s) read 
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(4.6e) 
(4.6f) 
(4.6g) 
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Using the above in Eq. (|2.10p . performing contraction of indices and truncating the resulting expression at the 2.5PN 
order, we can now write the total hereditary contribution at 2.5PN order, solely expressed as a function of our PN 
parameter 7 and it reads 
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(4.7) 



Note again, that leading order hereditary contribution (1.5PN tail) is proportional to various powers of s and hence 
would be absent when we specialize our result to case (b). The reason is similar to the one given at the end of 
Sec. IIV ATj to explain the absence of the Newtonian terms in instantaneous part for case (b). Observe that, the first 
two terms of Eq. (|2.10l) arc proportional to the I^"^^ and 1^*^^ , and these arc the ones which should contributing at the 

1.5PN order. But, since Newtonian order expression for vanishes for n > 2, for the case of infall from infinity, 
there would be no contribution at the 1.5PN order for case (b). 



2. Case (h): Infall from infinity 

(n) 

Using the argument, that the Newtonian order expression for II-^, vanishes for n > 2 in the case of infall from 
infinity, in Eq. (j2.10p , we can immediately see that only first two terms of Eq. (|2.10p are going to contribute to the 
linear momentum flux. And thus we need to evaluate only the integrals appearing in these two terms. In this case, 
the relevant integrals take the following form 
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With, required derivatives of the source multipole moments, expressed in terms of z, and the IPN trajectory (given 
by Eq. p.7p ) we can evaluate these integrals and they read 
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Using the above result in Eq. p.lOp . we can write the complete hereditary contribution at 2.5PN order, as a function 
of our PN parameter 7, and it reads 
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Total Linear Momentum Flux 



1. Case (a): Inf all from a finite distance 



For this case, Eq. (|4.ip and Eq. (|4.7p can be added to write the complete 2.5PN accurate expression for the linear 
momentum flux, expressed as a function of the parameter 7, and it reads 
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(4.11) 



2. Case (h): Inf all from infinity 



For this case, Eq. (|4.2p and Eq. (|4.10p can be added to get the complete 2.5PN accurate expression for the linear 
momentum flux, in harmonic coordinates, expressed as a function of the parameter 7, and it reads 
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V. RECOIL VELOCITY 

With, the 2.5PN expression for linear momentum flux emitted during the radial infall of two compact objects for 
two different situations (case (a) and case (b)), in harmonic coordinates, we can now use the momentum balance 
argument to write the loss rate of linear momentum from the source (through outgoing gravitational waves) and it 
reads 



du 



= -rp{u). 



(5.1) 
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The net loss of linear momentuni can be obtained by integrating the balance equation, i.e. 

AP' = - f duTpiu'). (5.2) 

J —oo 



A. Case (a): Infall from a finite distance 

In this case, Eq. (j5.2p can be written as 

ru(zf) 

AP' = - / duTpiu) 
dz 



i(z) 
Gm Pf 



J->(7). (5.3) 



as 7 = {Gm/c^z) and dz = —(Gmlc^^'^)d'y. Here, Zf denotes some final separation where we would like terminate 
our integral. Also, two limiting values of the parameter, 7, are 7^ ~ (Grajc^Zi) and 7/ = [Gm/ Zf). 

We can use the 2.5PN expressions for the linear momentum flux (Eq. (|4.1ip ') and for z (Eq. p.4p ') in the above 
integral to compute the total loss of linear momentum from the source during the radial infall from an initial separation 
of Zi (7i) to a final separation of Zf (7/). Since, linear momentum fiux given by Eq. (j4.1ip involves some integrals 
(Eq. (|4.6p ) which have to be computed numerically, we can not have an analytical expression for the total loss of the 
linear momentum from the source and thus need to be computed numerically. The corresponding recoil velocity can 
be computed as 

AV' = AFV?7z (5.4) 

where, m is the total mass of the system. We shall present our estimates for the recoil velocity for the case of infall 
from a finite distance in the next section where we shall discuss all our findings. 



B. Case (b): Infall from infinity 



In this case, the loss of linear momentum can be given by the integral 
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The 2.5PN expressions for the linear momentum fiux (Eq. (|4.12p ) and for z (Eq. p.Sp ) can be used in the above to 
compute the total loss in the linear momentum during the radial infall of the two objects for the case of infall from 
infinity. Next, Eq. (|5.4p can be used to compute the corresponding expression for the recoil velocity. We find for the 
2.5 PN recoil velocity, in harmonic coordinates, expressed in terms of 7 as 
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Y,=0.2; 2.5PN Model v=0.2; 2.5 PN Model 




FIG. 1: Recoil velocity as a function of the mass parameter v (left panel) and as a function of the post-Newtonian parameter 
7/ (right panel) has been plotted. The parameter, is known as symmetric mass ratio of the binary; the parameters 
7/ = {Gm/c? Zf) and 7; = {Gm/c^Zi) are the post-Newtonian parameters characterizing the final and initial separation 
of the two objects, respectively. For the plot in the left panel, the value of the parameter 7/ has been fixed to 0.2, which 
corresponds to the final separation of 5 Gm/c? between the two objects and then the recoil velocity as a function of the 
parameter v has been plotted. Similarly, for the right panel, the value of the parameter v has been fixed to 0.2 and recoil 
velocity as a function of the parameter 7/ has been shown. These plots (both in the left and the right panel) also compare 
recoil velocity estimates for four different situations based on the binary's initial separation: 7^=0.01, 0.02, 0.05, and 0.0 which 
correspond to the initial separation of the two objects of 100 Gm/c? , 50 Gm/c? , 20 Gm/c? , and 00 (infinite initial separation 
case), respectively. 

VI. DISCUSSIONS AND CONCLUSIONS 

The 2.5PN accurate expressions for the linear momentum flux emitted during the radial infall of two compact 
objects for two different situations (infall from some finite initial separation and infall from infinity), in harmonic 
coordinates, expressed in terms of the post- Newtonian parameter 7 (related to the separation of the two objects), has 
been given by Eq. (|4.1ip and Eq. (|4.12p . Next, we use these expressions to compute the associated recoil velocity of the 
source. Equation (|5.6p gives the 2.5PN accurate analytical formula for the recoil velocity accumulated till any epoch 
during the binary's evolution (within the validity of PN approximations), for the case of infall from infinity, and can be 
used to compute related numerical estimates for the recoil velocity. Since linear momentum flux formula (Eq. (j4.1ip ). 
for the case which assumes the infall from some finite initial separation, involves some integrals (Eq. (j4.6|) ) which can 
only be evaluated numerically, it is not possible to give analytical PN expressions for the accumulated recoil velocity 
for this case. Figures [T] and [2] show the numerical estimates for the recoil velocity accumulated during the radial infall 
of two compact objects and we shall discuss them one by one. 

Figure [T] plots recoil velocity as a function of v (left panel) and as a function of the parameter 7^ (right panel). 
Here, 7/ is our post-Newtonian parameter given by 7/ — (Gm/c? Zf). For the plots in the left panel of Fig. [T]the 
value of the parameter 7^ has been fixed to 0.2 and then the recoil velocity has been plotted as a function of v for 
the range of = 0.01 (nearly test particle limit) io v = 0.24 (nearly symmetric binary). The right panel shows the 
variations in recoil velocity estimates as a function of the parameter 7^ for a range of values between 7/ = 0.01 to 
7/ = 0.2, for a binary with v = 0.2. These plots (both in the right and the left panel) also compare the recoil velocity 
estimates for four different situations related to the initial separation of the two objects under the radial infall. The 
recoil velocity estimates have been plotted for four different values of the parameter 7^ = [Gm/c? Zi): 7^=0.01, 0.02, 
0.05 and 0.0 which correspond to the initial separation of 100 Gm/c? ^ 50 Gm/? ^ 20 Gm/? ^ and 00 (infinite initial 
separation case), respectively. 

Based on the estimates shown in the left panel of Fig. [TJ we find that the recoil velocity is maximum for a binary 
with V ^ 0.19 and is of the order of ^ 1.6kms~^. Also, the behavior of the plots is as one would expect: recoil 
velocity is maximum for the infinite initial separation case and estimates become smaller for situations which assume 
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FIG. 2: Recoil velocity as a function of the parameter v has been shown. For all the plots, the value of the parameter 7/ has 
been fixed to 0.2 (which corresponds to the final separation of 5 Gm/c? between the two objects under the radial infall). Plots 
in different panels also compare the results with different FN accuracy for four different situations: 7i=0.01, 0.02, 0.05, and 
0.0 which correspond to the initial separation (of the two objects in the problem) of 100 Gm/c? , 50 Gm/c? , 20 Gm/c? , and 
00 (infinite initial separation case), respectively. 



that infall shall proceed from smaller separations.^ However, we observe that estimates for the reeoil veloeity for all 
four situations (7,;=:0.01, 0.02, 0.05 and 0.0) are of the same order, indicating that most of the contribution comes 
from late stages of the infall. 

Although, we are not aware of a study which provides recoil velocity accumulated only during the premerger phase 
of a binary under the radial infall, a comparison with some other analytical/numerical work (which also involve 
contributions from the merger phase of the binary evolution) will be useful. For our purpose (head-on collision of 
two nonspinning compact objects), closest comparisons can be made using the results of [2^ (Numerical Relativity) 
and of [2g| (black hole perturbation theory). As compared to the recoil velocity estimates of about 2-5 kms~^ of [25| 
for a black hole binary (with v = 0.24) under radial infall, our estimates using (Eg . (j5.6p ) suggest a recoil velocity 
of the order of 0.95 kms~^ for the same system {i.e. with v = 0.24). Reference [2a] suggests that the recoil velocity 
accumulated during the head-on infall and plunge of a test particle in to a Schwarzschild black hole is given by 
IW = 8.73 X 10^^ V c, which, compared to our estimates of the recoil velocity using the test particle limit of Eq. (|5.6p 
(Al/ = 4.06 X 10~^!^c), is larger by a factor of two. The difference between our estimates and other related estimates 
is possibly due to the fact that we do not evolve our system till it merges. 

Figure [2] plots the recoil velocity as a function of v. For all the plots, the value of the parameter 7/ has been 
fixed to 0.2. Four panels correspond to the four initial separations which have been discussed above while describing 



Note that for finite separation cases (7^=0.01, 0.02, 0.05), initially the contribution exceeds as compared to the case of infinite initial 
separation (7^=0.0): this is not surprising since this contribution comes from the Newtonian terms which are absent in infinite initial 
separation case. 
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Figure [2l Each panel compares the recoil velocity estimates using results with different PN accuracy (Newtonian, 
••• , 2.5PN). It should be noted that we are terminating all our computations at = 0.2 [i.e when the distance 
between the two objects is 5 Gmjc?). The reason for this is related to the validity of our formulas beyond this final 
separation. Generally, it is believed that when higher order PN corrections start becoming comparable to the leading 
order contribution in the series and such a series becomes less reliable. A few checks with our analytical expressions 
indicate that these estimates are reliable for separations larger than 5 Gm/c^ (7 = 0.2) and this is why we terminate 
all our computations at this value (7 = 0.2). 

Appendix A: The 2.5PN linear momentum flux and recoil velocity in ADM coordinates 

In the above, we have given the 2.5PN accurate analytical expression for the linear momentum flux due to radial 
infall of two compact objects under mutual gravitational influence, in harmonic coordinates. In this section we shall 
provide equivalent formulas in ADM coordinates. 

1. Case (a): Infall from a finite distance 

The 2.5PN accurate analytical expression for the linear momentum flux in ADM coordinates can be obtained by 
using the following relation 

(-^p)adM = -^p + ^(Har-^ADM)-^^- (Al) 

Here, Tp is given by Eq. (|4.1ip and (5(Har^ADM)-7^p reads 

(A2) 

2. Case (b): Infall from infinity 

In this case, expression for the linear momentum flux in ADM coordinates can be obtained using Eq. (jAip . with 
T\, given by Eq. (|4.12p and (5(Har^ADM)-^p as 

'5(Har^ADM)-F> = JJ^t^^^ I ^ ^/Y^,,- ^ ^ ^^^^ 

In this case we can also write the recoil velocity expression in ADM coordinates by using the following relation 

(AF^)adM = A"^^'' + <5(Har^ADM)A1/\ (A4) 

Here, Ay is given by Eq. (|5.6|) and (5(Har-i.ADM) is given by 



-^(Har^ADM)^* - ^C7^Vl - ^Vv' ^ + 



1 . V 

6 



n,. (A5) 



Acknowledgments 

I thank Bala R. Iyer for suggesting this problem. I thank Bala R. Iyer and K. G. Arun for discussions and useful 
suggestions on the manuscript. 



[1] S. A. Hughes, M. Favata, and D. E. Holz, in Growing Black Holes: Accretion in a Cosmological Context , edited by A. Mer- 
loni, S. Nayakshin, & R. A. Sunyaev (2005), pp. 333-339, astro-ph/0408492. 
[2] D. Merritt, M. Milosavljevic, M. Favata, S. A. Hughes, and D. E. Holz, Astropliys. J. 607, L9 (2004), astro- ph/0402057. 



17 



[3] D. Richstone, E. A. Ajhar, R. Bender, G. Bower, A. Dressier, S. M. Faber, A. V. Filippenko, K. Gebhardt, R. Green, L. C. 

Ho, et al.. Nature 395, A14 (1998), astro-ph/9810378. 
[4] M. J. Fitchett, Mon. Not. Roy. Soc. 203, 1049 (1983). 
[5] A. G. Wiseman, Phys. Rev. D 46, 1517 (1992). 

[6] L. Blanchet, M. S. S. Qusailah, and C. M. Will, Astrophys. J 635, 508 (2005), astro-ph/0507692. 
[7] L. Kidder, Phys. Rev. D 52, 821 (1995). 

[8] E. Racine, A. Buonanno, and L. E. Kidder, Phys. Rev. D 80, 044010 (2009), arXiv:0812.4413 [gr-qc]. 

[9] M. Favata, S. A. Hughes, and D. E. Holz, Astrophys. J. 607, L5 (2004), astro-ph/0402056. 
[10] A. Buonanno and T. Damour, Phys. Rev. D 59, 084006 (1999), gr-qc/9811091. 
[11] T. Damour, Phys. Rev. D 64, 124013 (2001), gr-qc/0103018. 
[12] T. Damour and A. Gopakumar, Phys. Rev. D 73, 124006 (2006), gr-qc/0602117. 
[13] C. F. Sopuerta, N. Yunes, and P. Laguna, Astrophys. J. 656, L9 (2007), astro-ph/0611110. 

[14] C. F. Sopuerta, N. Yunes, and P. Laguna, Phys. Rev. D 74, 124010 (2006), Erratum-ibid.D75:069903,2007, Erratum- 

ibid.D78:049901,2008, astro-ph/0608600. 
[15] A. Le Tiec, L. Blanchet, and C. M. Will, Class. Quant. Grav. 27, 012001 (2010), 0910.4594. 
[16] C. K. Mishra, K. G. Arun, and B. R. Iyer, Phys. Rev. D 85, 044021 (2012), arXiv:1111.2701 [gr-qc]. 
[17] M. CampaneUi, Class. Quant. Grav. 22, S387 (2005), astro-ph/0411744. 
[18] J. G. Baker et al., Astrophys. J. 653, L93 (2006), astro- ph/0603204. 

[19] F. Herrmann, I. Hinder, D. Shoemaker, and P. Laguna, Classical and Quantum Gravity 24, S33 (2007). 

[20] J. A. Gonzalez, U. Sperhake, and B. Bruegmann, Phys. Rev. D 79, 124006 (2009), arXiv:0811.3952 [gr-qc]. 

[21] F. Herrmann, I. Hinder, D. Shoemaker, P. Laguna, and R. A. Matzner, Astrophys. J. 661, 430 (2007), gr-qc/0701143. 

[22] M. Koppitz et al., Phys. Rev. Lett. 99, 041102 (2007), gr-qc/0701163. 

[23] M. CampaneUi, C. O. Lousto, Y. Zlochower, and D. Merritt, Astrophys. J. 659, L5 (2007), gr-qc/0701164. 
[24] J. A. Gonzalez, M. D. Hannam, U. Sperhake, B. Bruegmann, and S. Husa, Phys. Rev. Lett. 98, 231101 (2007), gr- 
qc/0702052. 

[25] D.-I. Choi et al., Phys. Rev. D 76, 104026 (2007), gr-qc/0702016. 

[26] T. Nakamura and M. P. Haugan, Astrophys. J. 269, 292 (1983). 

[27] Z. Andrade and R. H. Price, Phys.Rev. D 56, 6336 (1997), gr-qc/9611022. 

[28] P. Anninos and S. Brandt, Phys. Rev. Lett. 81, 508 (1998), gr-qc/9806031. 

[29] D. A. Nichols and Y. Chen, Phys.Rev. D 85, 044035 (2012), arXiv: 1109.0081 [gr-qc]. 

[30] G. Lovelace, Y. Chen, M. Cohen, J. D. Kaplan, D. Keppel, et al., Phys.Rev. D 82, 064031 (2010), arXiv:0907.0869 [gr-qc]. 
[31] A. E. Simone, E. Poisson, and C. M. Will, Phys. Rev. D 52, 4481 (1995), gr-qc/9506080. 

[32] C. K. Mishra and B. R. Iyer, Phys. Rev. D 82, 104005 (2010), erratum-ibid. 84, 069908 (2011), arXiv: 1008.4009 [gr-qc]. 
[33] K. Thorne, Rev. Mod. Phys. 52, 299 (1980). 

[34] L. Blanchet, G. Faye, B. R. Iyer, and S. Sinha, Class. Quantum. Grav. 25, 165003 (2008), arXiv:0802.1249. 
[35] L. Blanchet and G. Faye, Phys. Rev. D 63, 062005 (2001), gr-qc/0007051. 
[36] L. Blanchet, Living Rev. Rel. 9, 4 (2006), gr-qc/0202016. 

[37] L. Blanchet, B. R. Iyer, and B. Joguet, Phys. Rev. D 65, 064005 (2002), Erratum-ibid 71, 129903(E) (2005), gr-qc/0105098. 
[38] L. Blanchet and B. R. Iyer, Phys. Rev. D 71, 024004 (2004), gr-qc/0409094. 

[39] L. Blanchet, T. Damour, and G. Esposito-Farese, Phys. Rev. D 69, 124007 (2004), gr-qc/0311052. 

[40] L. Blanchet, T. Damour, G. Esposito-Farese, and B. R. Iyer, Phys. Rev. D 71, 124004 (2005), gr-qc/0503044. 



